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Abstract 



We calculate the complete electroweak 0(a) corrections to p'p — > l + l~X(l = 
e, n) in the Standard Model (SM) of electroweak interactions. They comprise 
weak and photonic virtual one-loop corrections as well as real photon radia- 
tion to the parton-level processes qq — ► 7, Z — > l + l~ . We study in detail the 
effect of the radiative corrections on the l + l~ invariant mass distribution, the 
cross section in the Z boson resonance region, and on the forward-backward 
asymmetry, ^4fb , at the Fermilab Tevatron and the CERN Large Hadron Col- 
lider. The weak corrections are found to increase the Z boson cross section by 
about 1%, but have little effect on the forward-backward asymmetry in the 
Z peak region. Threshold effects of the W box diagrams lead to pronounced 
effects in ^4fb at m(l + l~) ~ 160 GeV which, however, will be difficult to ob- 
serve experimentally. At high di-lepton invariant masses, the non-factorizable 
weak corrections are found to become large. 
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I. INTRODUCTION 



Drell-Yan production in hadronic collisions, pp — ► l + l~X (I = e, //), is an interesting 
process for a number of reasons. Low mass Drell-Yan production is of interest because of the 
sensitivity to parton distribution functions (PDFs) at small x values In the Z-boson 

resonance region, measurement of the Z-boson mass, Mz, and width, Yz, and comparison 
with the values obtained at LEP helps to accurately calibrate detector components which 
is important for the determination of the W mass Measuring the forward-backward 

asymmetry, A^b, in the vicinity of the Z pole || makes it possible to extract the effective 
weak mixing angle. The ratio, R, of the W — > Iv and Z — > cross sections can be 
used to extract information on the width of the W boson |7U§. Finally, above the Z peak, 
one can search for physics beyond the SM such as extra neutral gauge bosons ||, effects of 
large extra dimensions [p!0| , p| , or composite quarks and leptons [12| in either the di-lepton 
invariant mass distribution or the forward-backward asymmetry. 

With the anticipated large data sets of 2 — 20 fb _1 for Run II of the Fermilab Tevatron and 
100 fb _1 per year at the CERN Large Hadron Collider (LHC), it is very important to fully 
understand and control higher order QCD and electroweak corrections to Drell-Yan produc- 
tion. A complete calculation of the full 0(a) radiative corrections to pp — > 7, Z — > l + l~ 
has not been carried out yet. In a previous calculation, only the QED corrections had been 
included [|H| while "genuine" weak corrections were ignored. At the Tevatron (LHC), the 



expected statistical uncertainty on the Z boson cross section for 2 fb _1 (10 fb _1 ) is approxi- 
mately 0.2% (0.05%) per lepton channel. In contrast, the genuine weak corrections in the Z 
peak region are expected to be of (9(1%) in magnitude and grow with the di-lepton invariant 
mass, m(l + l~), similar to di-fermion production in e + e _ collisions [T3]. It is thus necessary 



to include these corrections when data and SM prediction are compared. Furthermore, in 
order to properly calibrate the Z boson mass and width using the available LEP data, it 
is desirable to use exactly the same theoretical input which has been used to extract Mz 
and Yz at LEP, i.e. to include the reducible and irreducible O^mf/M^) corrections to 
the effective leptonic weak mixing parameter, sin 2 9 l eS , and the W mass, Mw [15] , in the 



calculation. Here, g is the SU(2) L coupling constant, and m t is the top quark mass. 

In this paper, we present a complete calculation of the electroweak 0{a) corrections to 
pp — > 7, Z — ► l + l~ which also takes into account the 0(g 4 m 2 /M^ v ) corrections to sin 2 ^ 
and Mw For the numerical evaluation, we use the Monte Carlo method for next-to- leading- 
order (NLO) calculations described in Ref. [T^]. With the Monte Carlo method, it is easy 
to calculate a variety of observables simultaneously and to simulate detector response. The 



QED corrections are taken from Ref. |13| and the collinear singularities associated with 



initial state photon radiation are removed by universal collinear counter terms generated 



by "renormalizing" the parton distribution functions fl7H19|l, in complete analogy to gluon 



emission in QCD. Final state charged lepton mass effects are included in our calculation in 
the following approximation. The lepton mass regularizes the collinear singularity associated 
with final state photon radiation. The associated mass singular logarithms of the form 
ln(s/m 2 ), where s is the squared parton center of mass energy and mi is the charged lepton 
mass, are included in our calculation, but the very small terms of 0(m 2 /s) are neglected. 

The technical details of our calculation are described in Sec. II. The electroweak cor- 
rections consist of the set of electroweak loop contributions, including virtual photons, and 
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of the emission of real photons. To regularize the ultraviolet divergences associated with 
the virtual corrections, we use dimensional regularization in the on-shell renormalization 
scheme ||20|| . After a brief summary of the calculation of the QED corrections [103], ana- 
lytical expressions for the genuine weak corrections are presented. In the Z pole region, 
the leading universal electroweak corrections can be expressed in terms of effective vector 
and axial vector couplings. These corrections can thus be taken into account in form of an 
effective Born approximation (EBA) where the tree level vector and axial vector coupling 
constants in the expression of the Born cross section are replaced by the effective vector and 
axial vector couplings. The remaining non-factorizable weak corrections are small in the Z 
pole region, but become important at high di-lepton invariant masses due to the presence 
of large Sudakov-like electroweak logarithms of the form \n(m(l + 1~) / My) {V = W, Z) ||14|| . 
In Sec. II, we also present a numerical comparison of the full 0{a) cross section and the 
forward-backward asymmetry at parton level with that obtained in the EBA and the Born 
approximation. Such a comparison is helpful to gain insight into how the weak corrections 
affect measurable quantities. 

Numerical results for pp collisions at y'i = 2 TeV and for pp collisions at y/s = 14 TeV 
are presented in Sec. III. When 0{g i ml/M^) corrections to sin 2 9 l eS and M w are taken into 
account, the weak corrections increase the Z boson cross section by about 1% but have little 
effect on Aps in the Z pole region. Threshold effects of the W box diagrams are found to lead 
to small but pronounced effects in the forward-backward asymmetry at m(l + l~) ~ 160 GeV. 
At large di-lepton invariant masses, the weak corrections reduce the differential cross section 
by 0(10%). The effect on Afb is somewhat smaller. Finally, our conclusions are presented 
in Sec. IV. 



II. ELECTROWEAK RADIATIVE CORRECTIONS TO NEUTRAL CURRENT 

DRELL-YAN PROCESSES 

The parton-level differential Born cross section for charged lepton pair production via 
photon and Z boson exchange in quark-antiquark annihilation (I = e,fi), 

q(p) + q{P) - 7, ^ - l + (k + ) + r(k-), (1) 

is given by 

dfr<°> = dP 2f ^EK + Azl'fr i, u) , (2) 

where the sum is taken over the spin and color degrees of freedom of the initial and final 
state fermions, and dPy denotes the two-particle phase space element. The factor 1/12 
results from summing and averaging over the quark spin and color degrees of freedom. The 
matrix elements A® and A° z describe the photon and Z boson exchange processes at lowest 
order in perturbation theory. In terms of the kinematical variables of the parton system 

s = (p + p) 2 , t=(p-k + ) 2 , u=(p-k„)\ (3) 

the squared Born matrix elements for massless external fermions are 
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1 2 

7 1 



2^2, ,2 (< 2 +« 2 ) 



8 (4vra) 

12 



£ 27£e(./4^*) = 16 g 9 g, a g a, (4vra) + w 2 ) - a y a ; (t 



7^ex(s) 



(4) 



with 



(i?-24Q/), 



r 3 



(5) 



parametrizing the Z// (/ = Z,g) couplings. Here, Qf and J| denote the charge and third 
component of the weak isospin quantum numbers of the fermion and s w = sin 9wi Cw = 
cos^vf with 6\y being the weak mixing angle, a = a(0) is the electromagnetic fine structure 
constant. The pole in the Z boson propagator is regularized by assuming a complex Z boson 
mass M r 



X(s) = Ana 



(s - Ml) ■ 



(6) 



In a perturbative calculation of the Z propagator, a Dyson resummation of one-particle- 
irreducible (1PI) (renormalized) Z self energies is performed. The imaginary part of M 2 is 
related to the Z decay width V z by unitarity while the real part is given by 



IZeMzis = Mi 



M 2 Z . 



(7) 



In the on-shell renormalization scheme |20| which we use, the physical Z boson mass Mz 
is equal to the renormalized mass. 

The Dyson resummation introduces the problem of defining the mass and width of the 
Z boson, and a gauge invariant description of the scattering amplitude, order-by-order in 
perturbation theory. As discussed in detail in Ref. J21j, the 0(a 2 ) contributions to the 
imaginary part of M 2 must be taken into account for a description of the Z resonance at 
the one-loop level. A consistent expansion of ImM c 2 to 0(a 2 ) yields 



ImM 2 (s) 



l+TZeil z {M 2 z) +t 



J (2) 



.S - 



(8) 



£ z and E^-j in Eq. d§|) are the transverse parts of the renormalized one and two-loop cor- 
rected 1PI Z self energies. fl z (see Eq. fl]JD) and £ z >^ z (see Eqs. flBl]) - H)) are the 
(renormalized) self energy insertions into the Z and photon propagators. The last term in 
Eq. (U) takes into account that the photon and Z boson do not propagate independently 
beyond leading order in perturbation theory. The evaluation of ImM 2 at s — M§ corre- 
sponds to a Laurent expansion around the complex pole and leads to a description of the Z 



resonance in terms of a constant width, i.e. M 2 



M 2 



iM z T 



zy z 



(0+1) 



At LEP energies the 



s-dependence of the one and two loop contributions to the imaginary part of the 1PI Z self 
energy can be approximated by 
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lm± z (s) = Tm± z {M 2 z ) Xmtf 2) (s) = lm£f 2) (M 2 z ) (9) 
z z 

so that M 2 = M§ — i(s / Mz)T^ +1 \ Both descriptions are related by a transformation of the 
parameters of the Z resonance, Mz and Tz, and the residue of the complex pole, and are 
thus equivalent p2| . In the following, we use the s-dependent width approach. The one-loop 
corrected Z boson decay width, , is discussed in Appendix |A]. 

The electroweak 0{a) corrections to neutral-current Drell-Yan processes naturally de- 
compose into QED and weak contributions, i.e. they form gauge invariant subsets, and thus 
can be discussed separately. The observable NLO cross section is obtained by convoluting 
the parton cross section with the quark distribution functions f q /A(%,Q 2 ) {s = Xix 2 s) and 
summing over all quark flavors q 

d(T ( s ) ='E£dx 1 dx 2 (f q/A (xi,Q 2 )fg/B(x 2 ,Q 2 ) + (q^ qj\ [d&(° +1 \s,t,u) (10) 



u 



with (A,B) = (p,p) for the Tevatron and (p,p) for the LHC. da^ 0+1 ^ comprises the NLO 
cross section including the weak corrections and g^ctqed describes the QED part, i.e. virtual 
and real photon emission off the quarks and charged leptons. The parton distribution 
functions depend on the QCD renormalization and factorization scales, \i r and fif, which 
we choose to be equal, \i r = p,f = Q. The radiation of collinear photons off quarks requires 
the factorization of the arising mass singularities into the PDFs which introduces a QED 
factorization scale /iqed as will be explained in more detail in the next section. 



A. QED corrections 

QED radiative corrections consist of the emission of real and virtual photons off the 
quarks and charged leptons. The 0(a) QED corrections to qq — >• 7, Z — ► l + l~ can be 
further divided into gauge invariant subsets corresponding to initial and final-state radiation. 
Since the incoming quarks are assumed to be massless in the parton model, initial state 
photon radiation results in collinear singularities. These singularities are universal to all 
orders in perturbation theory and can be absorbed by a redefinition {renormalization) of 
the PDFs |17| . This can be done in complete analogy to the calculation of QCD radiative 
corrections. As a result, the renormalized parton distribution functions become dependent 
on the QED factorization scale /xqed which is controlled by the well-known Gribov-Lipatov- 



Altarelli-Parisi (GLAP) equations f23|. These universal photonic corrections can be taken 



into account by a straightforward modification [18,19] of the standard GLAP equations 
which describe gluonic corrections only. The modification consists of an additional term 
which is proportional to the electromagnetic fine-structure constant, a, resulting in modified 
distribution functions qf(x, /Uq ED ) for quarks with flavor /. The gluon distribution g(x, /Uq ED ) 
is only affected indirectly by QED corrections through terms of 0(aa s ). The QED induced 
terms in the GLAP equations lead to small, negative corrections at the per-mille level to the 
distribution functions for most values of x and /^qed [ 24 ] . Only at large x ~ 0.5 and large 



yUq ED ~ 10 GeV do the corrections reach the magnitude of one per cent. 
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In order to treat the 0(a) initial state QED corrections to Drell-Yan production in 
hadronic collisions in a consistent way, QED corrections should be incorporated in the global 
fitting of the PDFs, i.e. all data which are used to fit the parton distribution functions should 
be corrected for QED effects. Current fits ||25lj26| to the PDFs do not include QED correc- 
tions. The missing QED corrections introduce an uncertainty which, however, is probably 
much smaller than the present experimental uncertainties on the parton distribution func- 
tions. 

Absorbing the collinear singularity into the PDFs introduces a QED factorization scheme 
dependence. The squared matrix elements for different QED factorization schemes differ by 
the finite 0(a) terms which are absorbed into the PDFs in addition to the singular terms. 

For our calculation, we have taken the 0(a) QED corrections from Ref. ||13|| . The 



calculation presented in Ref. is based on an explicit diagrammatic approach. The 
collinear singularities associated with initial-state photon radiation are factorized into the 
parton distribution functions as described above, but QED corrections are not taken into 
account into the GLAP evolution of the PDFs. The 0(a) QED corrections are implemented 
both in the QED MS and DIS schemes, which are defined analogously to the usual MS p7[ 
and DIS p|[ schemes used in QCD calculations. All numerical calculations in this paper 
are performed using the QED DIS scheme. The QED DIS scheme is defined by requiring 
the same expression for the leading and next-to-leading order structure function F 2 in deep 
inelastic scattering. Since F 2 data are an important ingredient in extracting PDFs, the effect 
of the 0(a) QED corrections on the PDFs should be reduced in the QED DIS scheme. The 
collinear singularities associated with photon radiation from the final state lepton lines are 
regulated by finite lepton masses. 



B. Non-QED corrections and the effective Born approximation 

The non-QED corrections consist of self-energy contributions to the photon and Z prop- 
agators, vertex corrections to the r y/Z-l + l~ and j/Z-qq couplings (see Fig. |l]), and box 
diagrams with two massive gauge bosons (see Fig. 0). Since we neglect all non-logarithmic 
fermion mass terms, there is no Higgs boson contribution to the box diagrams and vertex cor- 
rections. The calculation of the radiative corrections is performed in the 't Hooft-Feynman 
gauge. To regularize and remove the arising ultraviolet (UV) divergences we use dimensional 



regularization in the on-shell renormalization scheme as described in Ref. [20 



In the following we closely follow Refs. PU[ and (RUf, in particular for the treatment of 



additional higher-order corrections. For an accurate description of the Z resonance, it is 
important to take these corrections into account. The NLO differential cross section at the 
parton level including weak 0(a) and leading 0(a 2 ) corrections is of the form 

da^ = dP v 1 £ |4 0+1 ) + 4 0+1) | 2 (s,t,u) + da hox (s,i,u) • (11) 

Here, dow describes the contribution of the box diagrams shown in Fig. |[ The matrix 
elements A^ 1 ^ comprise the Born matrix elements, A® z , the 7, Z and 7Z self energy 
insertions, including a leading log resummation of the terms involving the light fermions, 
and the one-loop vertex corrections, as shown in Fig. |l|. A^z can be expressed in terms of 
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effective vector and axial- vector couplings, gy'^'^ (f = l,q), such that the squared matrix 
elements for massless external fermions can be expressed in the form 



E I4 0+ 



1)12 



(47T«) 2 



X 



[l + neW{s)] 2 P 

8[(l^r + bl'?)(l^'T + bl'T) (P + u 2 ) 

AKe{g v \g^)Ke{g v «{gyy) (£ 2 - « 2 )1 , 



£1^ 



(0+1) 1 2 

z 



x 



[i + ^eiF(s)] 2 s 2 
8[(I^I 2 + I^I 2 )(I^T + I^'T) (P + e) 

±ne{g*\g z /y) TZe^g^T) (P - u 2 ) , 



^2^e(4 0+1) 4 0+1 )*) 



(4tt«) 



16Ke (x _1 (s) x 



[1 + KeW(s)} [1 + ftelF (s)] s 2 

(gl'W)* + &\g\*Y)(&"{gVY + 9v q teXTX* 2 - u 2 ) 



with 



9 Z /(S) =a f + G z /(s) , 
g Z y J {s)=v f + F^{s)+Q f 



W z (s) 

1 + W(s) ' 



,7,// 



<&{§) =Qf- F^(s) 



v 



and G 



A 



(12) 



(13) 



(14) 



(15) 



denote the renormalized vector and axial-vector formfactors which 



parametrize the weak corrections to the (j,Z)ff vertices. Il x ,X = 7, Z, jZ, describe 
the renormalized photon, Z and (7, Z) self energy insertions, 



TP(s) = 

rF z (s) = 

fl Z (s) = 



S 7 (g) 
s 

±~* z (s) 
s 

1 

s-M 2 z 



(^ Z (§)) 2 ' 

s + Et(s). 



(16) 
(17) 
(18) 



The box contribution da^ ox cannot be absorbed in effective couplings. However, in the Z 
resonance region, the box diagrams can be neglected and the NLO cross section da^ 0+1 ^ is of 
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Born-structure. For more details about the self energies, form factors and box contributions 
we refer to Appendices [B| and 0, respectively. In Appendix [B] we also describe the inclusion 
of leading higher-order (irreducible) QCD and electroweak corrections connected to the p 
parameter. 

In the Z resonance region, the dominant contributions to the non-photonic electroweak 
corrections can be taken into account by redefining several quantities appearing in the expres- 
sion of the Born cross section given in Eq. (0). In the resulting effective Born approximation 
(EBA), the fine structure constant, a, is replaced by the running electromagnetic coupling, 

OL 

a -> a(s) = - — — , Aor(s) = -1leTli rm (s) , (19) 

1 — Aa{s) 

where ffj erm denotes the fermion-loop contribution to the photon vacuum polarization. x(s) 
is expressed in terms of the physical W and Z boson masses, the Z width measured at LEP, 
and the Fermi constant G^, 

X(*) = W5S,^ + * ^ , (20) 



where 



s 



, 2 = / - _ Mw 
w I M 2 



1-lfJ- (21) 



Finally, the vertex and self energy corrections are taken into account by replacing 

v f - vf = -L ^ {If - 2Q f sin 2 0f s ) , / = I, q , (22) 



IV 



where 



^^h-pmi), (23, 



sin 



410/1 V Keg z /(MI 



is the effective electroweak mixing parameter for leptons at the Z peak, sin 9 eS has been 
measured at LEP and the SLC The effective weak mixing parameter for up and down type 
quarks is approximately given by 

sin 2 9l s w sin 2 6 l cS - 0.0001, (24) 
sin 2 d cS w sin 2 6 l cS - 0.0002. (25) 

Above the Z peak region, the effective Born approximation becomes insufficient for two 
reasons: the effective couplings are not static but grow as functions of s, and the box 
diagrams are no longer negligible. Their contribution increases strongly with energy and 
they contribute significantly at high invariant masses of the lepton pair. 
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C. Numerical discussion of the non-QED corrections at the parton level 



Before we present results at the hadron level, we discuss the observables of interest at 
the parton level, where many characteristics of the weak corrections manifest themselves. 
For the numerical evaluation we chose the following set of SM input parameters [j3l[] : 



1.16639 x 1(T 5 GeV~ 2 , 


a = 


= 1/137.0359895, 


M z = 91.1867 GeV, 


a s 


= a s {M 2 z ) = 0.119, 


m e = 0.51099907 MeV, 




= 0.105658389 GeV, m T = 1.777 GeV, 


m u = 0.0464 GeV, 


m c 


= 1.5 GeV, m t = 174 GeV, 


m d = 0.0465 GeV, 


m s 


= 0.15 GeV, m fe = 4.7GeV. 



(26) 

The fermion masses only enter through loop contributions to the vector boson self energies 
and as regulators of the collinear singularities which arise in the calculation of the QED con- 
tribution. Non-zero light quark masses are only used in the calculation of the vector boson 
self energies. The light quark masses are chosen such that the value for the hadronic contri- 
bution to the photon vacuum polarization for five active flavors, Aa^ d (M|) = 0.028 [32 



which is derived from low-energy e + e~ data with the help of dispersion relations, is recovered. 

The W mass and the Higgs boson mass, Mh, are related via loop corrections. A 
parametrization of the W mass which, for 65 GeV < Mu < 1 TeV, deviates by at most 
0.4 MeV from the theoretical value including the full fermionic two-loop contributions is 



given in Ref . |33[ . Here we use the parametrization of Ref . pi 



M W = M° V - 0.0581 In ( —^—) - 0.0078 In 2 (—M^—) 



0.085 



at. 



V100 GeV/ ' V100 GeV/ ' V0.118 

_ aiU (^_ 1 ) +aM7(( _^_ ) ._ l) (27) 

with Myy = 80.3805 GeV, which was used in the analysis of the LEP data. The parametriza- 
tion of Eq. (fZ7|) reproduces the result of Ref. to 0.2 MeV for 75 GeV < M H < 350 GeV. 
For the numerical discussion we choose 

Mh = 120 GeV, (28) 

which is consistent with current direct f35| and indirect bounds |36[[ , and work in the s- 



dependent width scheme. The Z-boson decay width is calculated including electroweak 
and QCD corrections as described in Appendix [A|. The NLO prediction for the Z boson 
width is also used in the calculation of the lowest-order and EBA predictions. For the 
input parameters listed in Eq. ( ^6|) we obtain for the effective leptonic weak mixing angle of 
Eq. (H) sin 2 9 l cS = 0.23167, and T z = r§ +1) = 2.4932 GeV for the width of the Z boson. 



In the following we discuss the impact of the weak corrections on the total cross sections 
and the forward-backward asymmetries for uu (dd) — > 7, Z — > e + e~ as a function of the 
parton center of mass energy Almost identical results are obtained for the final 
state. We compare the full NLO result, dcr^ 0+1 ^ of Eq. ([11]), with the Born prediction, ci<T (0) 
of Eq. (0), and the result obtained in the EBA. In Fig. |3|, we show the relative corrections 
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a^ 0+1 >/ao — 1, in per cent, to the total cross sections for the processes uu — > e + e~ and 
dd — > e + e~ in the Z resonance region (Fig. |3]a), and at high parton center of mass energies 
(Fig. |3|b). o"o is either taken to be the Born cross section, a^°\ or the EBA prediction for the 
total cross section, oeba- The weak corrections are seen to enhance the total cross section 
by 10 — 12% below, and by 5 — 7% at and above the Z peak. The kink at y/§ ~ 350 GeV 
in the solid and dotted lines is due to the top quark threshold in the running coupling, 
a(s). Since this enhancement is mainly caused by universal electroweak corrections, i.e. 
the running of a and corrections connected to the p parameter, the EBA represents a good 
description of the NLO result in the Z resonance region (a^ 0+1 ^ /a — 1 < 1%). The difference 
between the full NLO result and that obtained in the EBA is a measure for the effects of 
the non-universal corrections. At higher parton center of mass energies the strong deviation 
of the EBA cross section from the NLO result is essentially due to the contribution of the 
box diagrams and to a lesser extent due to the energy dependence of the effective couplings. 
The deviation of the EBA cross section from the NLO result at large values of is the 
result of large Sudakov-like electroweak logarithms of the form \n(s/My) (V = W, Z) [|14|| . 
It is much more pronounced for uu — > e + e~ than for dd — > e + e~ which is due to the different 
Feynman diagrams which contribute to the two initial states. In addition to the box diagrams 
with Z boson exchange, only the crossed (direct) W box diagram contributes for up-type 
(down-type) quarks in the initial state (see Appendix Q). As a result, for dd — > e + e~ , the 
contribution from the box diagrams is less significant and the deviation of the EBA cross 
section from the NLO result is not as large. 

The forward-backward asymmetry at the parton level is given by 

a F (s) - a B (s) ~ f ( r) / ~\ r^-W^da,^ 

A FB = , p K , B )J , <7 (s) = / dt , 29 

o b (s) + a B (s) J-s/2(-s) dt 

with da = da^°\ d&EBA, da^ 1 ^ yielding the Born, EBA and NLO predictions of Afb, re- 
spectively. In Fig. f| we show the differences A^ — Ap^ n and Ap^° — Ap^ A for the parton 
level processes uu — > e + e~ and dd — > e + e~ in the vicinity of the Z-peak (Fig. f|a) and at 
large parton center of mass energies (Fig. |b). In the Z peak region, Ap-g is strongly reduced 
by the weak corrections. The reduction is more pronounced for the dd subprocess. In the 
vicinity of the Z boson peak, the EBA again provides a very good approximation. Due to 
the s-dependence of the effective couplings and the contributions from the box diagrams, 
the difference between Apg° and v4pi A rapidly increases in magnitude for \/§ > 100 GeV. 
At very high values of y/§, the weak corrections to uu — > e + e~ (dd — > e + e~) considerably 
diminish (enhance) the forward-backward asymmetry. 

As can been seen from Fig. ||b, the weak corrections also lead to a sharp peak (dip) in 
Afb for uu — > e + e~ (dd e + e~) in the vicinity of the W + W~ threshold, ~ 160 GeV. 
The peak/dip is due to the V — A nature of the coupling of the W to fermions and a 
threshold effect in the WW box and vertex diagrams. Figure |^ shows 5A FB = Ap^° — 
Ap^ A in more detail in the region around the W + W~ threshold. The solid lines show 
5A-pB including the full set of Feynman diagrams contributing to the non-photonic weak 
corrections. Disregarding the WW box diagrams removes a large portion of the peak/dip 
(dashed lines). The remaining effect is due to initial and final state vertex corrections 
involving two virtual W bosons. For uu — ► e + e _ these interfere destructively whereas for 
dd —>■ e + e~ there is constructive interference. The effect of the initial and final state vertex 
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corrections involving two virtual W bosons therefore is more pronounced for d-type quarks 
in the initial state. Unlike the W box diagrams, the Z box graphs have only a very small 
effect on the forward-backward asymmetry. This is due to the small vector coupling of the 
Z boson to the charged leptons. The effect of removing the Z box diagrams from the set 
of Feynman diagrams describing the non-photonic weak corrections is shown by the dotted 
lines in Fig. |5|. 



III. PHENOMENOLOGICAL RESULTS 
A. Preliminaries 

We shall now discuss the phenomenological implications of the 0(a) genuine weak cor- 
rections to di-lepton production at the Tevatron (pp collisions at = 2 TeV) and the 
LHC (pp collisions at yfs = 14 TeV). We first discuss the impact of the non- universal weak 
corrections on the lepton pair invariant mass distribution and the total cross section in the 
Z pole region. We then consider how the forward-backward asymmetry, A-ps, is affected 
by these corrections. The universal weak corrections are taken into account in form of the 
effective Born approximation described in detail in Sec. [II B| . The SM parameters used in 



our numerical simulations are listed in Eqs. fl2"6|) - fl2"8p . To compute the cross section we use 



the MRSR2 set of parton distribution functions ||37|| ) and take the renormalization scale, fi, 
and the QED and QCD factorization scales, /zqed and /iqcd, to be /i 2 = //q ED = /Uqcd = 

To simulate detector acceptance, we impose the following transverse momentum (p T ) 
and pseudo-rapidity (rf) cuts (7 = e, n) 

Pr(0>20GeV, \rj(l)\ < 2.5. (30) 



These cuts approximately model the acceptance of the CDF II [|38| and D0 [|39] detectors at 
the Tevatron, and the ATLAS and CMS ||l|] detectors at the LHC. Uncertainties in the 



energy measurements of the charged leptons in the detector are simulated in the calculation 
by Gaussian smearing of the particle four-momentum vector with standard deviation a which 
depends on the particle type and the detector. The numerical results presented here were 
calculated using a values based on the CDF II and ATLAS specifications. 

The granularity of the detectors and the size of the electromagnetic showers in the 
calorimeter make it difficult to discriminate between electrons and photons with a small 
opening angle. In such cases we recombine the four-momentum vectors of the electron and 
photon to an effective electron four-momentum vector. The exact recombination procedure is 
detector dependent. For calculations performed at Tevatron energies we recombine the four- 
momentum vectors of the electron and photon to an effective electron four-momentum vector 
if both traverse the same calorimeter cell, assuming a calorimeter segmentation of Ai] x A<p = 
0.1 x 15° (<fi is the azimuthal angle in the transverse plane). This procedure is similar to that 
used by the CDF Collaboration in Run I. The segmentation chosen corresponds to that of 
the central part of the Run I CDF calorimeter. At LHC energies, we recombine the electron 
and photon four-momentum vectors if their separation in the pseudorapidity - azimuthal 
angle plane, 



AR(e, 7 ) = v /(A^(e, 7 )) 2 + (A0(e, 7 )) 2 , (31) 
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IS 

Ai?(e,7) < R c = 0.07, (32) 
similar to the resolution expected for ATLAS [[|(J. Recombining the electron and photon 



four-momentum vectors for small opening angles of the two particles greatly reduces the 



effect of the mass singular logarithmic terms associated with final state photon radiation |13 

Muons are identified in a hadron collider detector by hits in the muon chambers. In 
addition, one requires that the associated track is consistent with a minimum ionizing par- 
ticle. This limits the energy of a photon to be smaller than a critical value E2 for small 
muon - photon opening angles. At the Tevatron we impose a E 1 < E] = 2 GeV cut for 



photons traversing the same calorimeter cell as the muon. At the LHC, following Ref. [40 



we require the photon energy to be smaller than E] = 5 GeV if Ai?(/i, 7) < 0.3. The cut on 



the photon energy increases the size of the QED corrections for m(n + /j,~) > 100 GeV [13 

We impose the cuts and lepton identification requirements described above in all subse 
quent numerical simulations, unless explicitly noted otherwise. 



B. Weak corrections to the di-lepton invariant mass distribution and the Z boson 

cross section 

QED corrections are known to have a profound impact on the shape of the di-lepton 



invariant mass distribution ||13|| . Due to the mass singular terms associated with final state 
photon radiation, the differential cross section is reduced in the Z peak region by about 
10% for electrons, and by about 20% for muons in the final state. Below the Z resonance 
region, final state photon radiation enhances the cross section by up to a factor 1.5 with the 
maximum effect occurring at m(/ + /~) 75 GeV. For m(l + l~) > 100 GeV, QED corrections 
reduce the e + e~ (fi + fi~) differential cross section by about 5% (12 — 15%). In contrast to final 
state photon radiation which significantly changes the shape of the di-lepton invariant mass 
distribution, initial state QED corrections are uniform and small (~ +0.4%). The distortion 
of the Breit-Wigner shape of the Z resonance curve due to final state QED corrections causes 
the Z boson mass extracted from data to be shifted by about —100 MeV for electrons, and 
—300 MeV for muons, in the final state |3]|| . 

Based on the results obtained at the parton level (see Sec. |U G| ), one expects that the non- 
universal 0(a) weak corrections are small in the vicinity of the Z resonance. Figure |6| shows 
the ratio [da°^ a ' / 'dm(/ + Z _ )]/[<i<7^ ED / 1 dm(l + l~)] at the Tevatron, where a 0<y01 ' denotes the 
full NLO cross section, and cr^ ED represents the cross section which includes the factorizable 
electroweak corrections in form of the effective Born approximation together with the 0(a) 
QED corrections. Very similar results are obtained at the LHC. For m(/ + Z~) < 50 GeV, the 
non-universal corrections are very small and negative. In the Z peak region, they enhance 
the differential cross section by up 1.2%. Finally, for m{l + l~) > 130 GeV, the non-universal 
weak corrections become negative and rapidly increase in magnitude. The small differences 
between the results for electrons and muons in the final state are mostly due to the different 
lepton identification requirements. The slight dip visible at m(l + l~) ~ 160 GeV is caused 
by the W pair threshold effect discussed in Sec. |11 Q 

Since the non-universal weak corrections are not uniform in the Z peak region, they are 
expected to shift the Z boson mass extracted from data upward by several MeV. This is 
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much smaller than the effect caused by QED corrections but may not be negligible in future 
hadron collider experiments. In order to make a more quantitative prediction of the shift in 
Mz due to the weak corrections, detailed simulations which fully take into account detector 
response need to be performed. These are beyond the scope of this paper. 

The ratio, R, of the W —>■ Iv and Z — > cross sections can be used to extract informa- 
tion on the width of the W boson Since the QCD corrections to W and Z production 
are very similar, they cancel almost perfectly in the W to Z cross section ratio; the 0(a s ) 
corrections to R are of (9(1%) or less, depending on the set of parton distribution functions 
used [|42|| . In addition many experimental uncertainties such as the luminosity uncertainty, 
cancel in the cross section ratio. Accurate knowledge of how electroweak corrections af- 
fect the W — > lu and the di-lepton cross sections in the Z resonance region is thus very 
important. 

The effect of the non-universal weak corrections on the invariant mass distribution 
is also reflected in the total cross section in the Z resonance region. In Table | we list the 
cross section ratios 

= . ( 33 ) 

("EW K-factor") and 

^QED 

("QED ^-factor") for 75 GeV < m(l+r) < 105 GeV (/ = e, \i) at the Tevatron with 
and without taking the cuts and lepton identification requirements described in Sec. [Ill A 



into account. Similar results are obtained at the LHC. One observes that the genuine 
weak interactions increase the cross section by about 1.0%. Approximately one-half of 
the enhancement is due to the (9(g 4 m 2 /M^) corrections to sin 2 9 l eS and M w . In contrast, 
the QED corrections decrease the cross section. The size of the QED corrections to the 
cross section depends on the flavor of the final state lepton and whether cuts and lepton 
identification requirements are taken into account or not. Without detector effects taken 
into account, the QED corrections are numerically more important than the genuine weak 
corrections. Due to the mass singular terms associated with final state photon radiation, 
the QED corrections for the e + e~ final state are larger than those in the muon case. The 
full electroweak corrections reduce the e + e~(7) (fi + jJT (7)) cross section in the Z resonance 
region by about 5% (2%). 

The recombination of electron and photon momenta when the opening angle between the 
two particles is small strongly reduces the effect of the QED corrections to the integrated 
e + e~ cross section. As a result, the effects of the QED corrections and the genuine weak 
corrections partially cancel. The net effect of the electroweak corrections is a decrease of 
the cross section by 1.2%. In the muon case, lepton identification requirements increase the 
magnitude of the QED corrections, and the full electroweak corrections decrease the cross 
section by more than 6%. 

The effect of the non-factorizable weak corrections on the di-lepton invariant mass distri- 
bution at the Tevatron for large values of m{l + l~) is shown in Fig. [7] where we plot the ratio 
of the complete 0(a 3 ) electroweak and the EBA differential cross section as a function of 
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m{l + l~) (dashed lines). In order to make the effect of the non-factorizable weak corrections 
more transparent, we also show the corresponding ratio for the case where only the O{o) 
QED corrections and the factorizable electroweak corrections in form of the effective Born 
approximation are taken into account (solid lines). Figure || shows the corresponding results 
for the LHC Due to the recombination of electrons and photons, the QED corrections reduce 
the e + e~ differential cross section by only 3 — 5% over the invariant mass regions considered. 
In the muon case, the cut on the photon energy for photons which have a small opening 
angle with the muon reduces the hard photon part of the 0(a 3 ) fi + fi~( , y) cross section. 
As a result, the QED corrections are much more pronounced and display a much stronger 
dependence on the di-lepton invariant mass than in the e + e~ case. The non-factorizable 
weak corrections are seen to increase rapidly in size with m(l + F). As mentioned before (see 
Sec. [11 C| ), this is due to the presence of Sudakov-like electroweak logarithms of the form 
\n{m{l + F) / M v ) (V = W, Z). Most of the effect is caused by the contribution of up-type 
quarks in the initial state (see Fig. |^). For m(e + e _ ) = 500 GeV (m(e + e~) = 1.5 TeV), 
the electroweak corrections reduce the cross section by about 10% (15%) at the Tevatron 
(LHC). In the muon channel, the 0(a) electroweak corrections are larger in magnitude 
than the 0(a s ) QCD corrections for /i + yU~ invariant masses larger than about 500 GeV. 
For m(fi + fi~) = 2 TeV at the LHC, the 0(a) electroweak radiative corrections reduce the 
cross section by more than 35%, which is approximately equal to the expected statistical 
uncertainty in a 200 GeV bin centered at m(fi + fi~) = 2 TeV for 100 fb -1 . It will thus be 
important to take into account the non-factorizable weak corrections when measuring the 
Drell-Yan cross section at large di-lepton invariant masses at the LHC. 

The results shown in Fig. |8| should be interpreted with caution. Since the non-factorizable 
weak corrections become large for di-lepton invariant masses above 1 TeV, they need to be 
resummed in order to obtain accurate predictions in this phase space region (for a recent re- 
view of the resummation of electroweak Sudakov-like logarithms see Ref. ||43||). A calculation 
of Drell-Yan production in hadronic collisions which includes resummation of electroweak 
logarithms has not been carried out yet. 



C. Weak corrections to the forward-backward asymmetry 



We now discuss how the non-universal weak corrections affect the forward-backward 
asymmetry, A^b- The expressions used in the literature to define A^b at the Tevatron |44j 



and the LHC |13|] are slightly different. For pp collisions at Tevatron energies, A^b usually 
is defined by 

F- B 



A 



FB 



where 



— dcos6*, 

o a cos 9* 



F + B 



B 



o da 



i dcos9* 



dcos8*. 



Here, cos 6* is given by [f44| , f45" 

COS0* 



m(l+F)^Jm 2 (l+F) + p^(l+F) 



P + (r)p-(i + )-p-(r)p + (i + ) 



(35) 



(36) 



(37) 
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with 



P ± = ^(E± Pz ), (38) 

where E is the energy and p z is the longitudinal component of the momentum vector. In this 
definition of cos6*, the polar axis is taken to be the bisector of the proton beam momentum 
and the negative of the anti-proton beam momentum when they are boosted into the 
rest frame. In pp collisions at Tevatron energies, the flight direction of the incoming quark 
coincides with the proton beam direction for a large fraction of the events. The definition 
of cos#* in Eq. (|3T| ) has the advantage of minimizing the effects of the QCD corrections 
(see below). In the limit of vanishing di-lepton p?, 0* coincides with the angle between the 
lepton and the incoming proton in the rest frame. 

QED corrections are known to have a significant effect on the forward backward asym- 
metry for 50 GeV < m(l + l~) < 90 GeV but are small for di-lepton masses larger than 
100 GeV 0. The difference 

AA FB = A FB (full EWK) - A F b(QED) (39) 

is the quantity which best displays how the non-universal weak interactions influence the 
forward-backward asymmetry. Here, ApB(full EWK) is the forward-backward asymmetry 
calculated taking the full 0(a) electroweak corrections and the 0(g 4 mf/My V ) corrections 
into account. v4fb(QED)j on the other hand, only includes the 0(a) QED corrections, 
in addition to the factorizable corrections absorbed in the EBA. Fig. |9| shows AAfb for 
di-lepton masses between 40 GeV and 200 GeV at the Tevatron. It demonstrates that the 
weak corrections have only a small effect on the forward-backward asymmetry for m(/ + /~) < 
200 GeV. The peak in Ay4 FB located at m{l + l~) « 160 GeV originates from threshold effects 
associated with the W box diagrams (see Figs. |4]b and |5[). Although the contributions from 
up- and down-type quarks in the initial state tend to cancel, a significant effect remains since 
the uu parton luminosity is much larger than that for dd pairs for the di-lepton invariant 
mass range of interest. The slight differences between electron and muon final states in 
AAfb originate from the different detector resolutions for the two final states. The peak 
in AAfb located at m{l + l~) « 80 GeV is also associated with the WW box diagrams. 
In this region one of the W bosons in the loop is on-shell, causing a small resonance like 
enhancement in the forward-backward asymmetry. 

The peak located at m{l + l~) « 160 GeV is a characteristic signature of the non- 
factorizable weak interactions in Drell-Yan production and it is interesting to investigate 
whether it may be observable in Run II. Since the size of the effect is small, one has to 
worry about how higher QCD corrections, detector effects and background processes af- 
fect the peak. QCD corrections are uniform in the region of interest and do not modify 
the structure of the peak [|6[]. Detector resolution effects broaden the peak and reduce its 



height. These effects are taken into account in Fig. || Higher order Coulomb corrections are 
also expected to modify the shape of the peak. Finally, backgrounds from ti — > l + l~vi?bb, 
W + W~ — ► l + l^vv and ZZ — > l + l~vi) have to be taken into account. The it and W + W~ 
backgrounds can either be subtracted using the experimentally determined efi^> T + X cross 
section, or suppressed by imposing missing transverse momentum and jet veto cuts. Re- 
quiring that no jets with Pt(J) > 20 GeV and \f](j)\ < 3.5 are observed and imposing a 
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f T < 20 GeV cut suppresses the ZZ and it backgrounds to negligible levels at Tevatron en- 
ergies. AAfb for the e + e~ final state including the contribution of the W + W~ background 
is shown by the dotted line in Fig. |^. The W + W~ background is completely negligible for 
m(e + e~) < 100 GeV. For larger invariant masses it slightly decreases AA FB but does not 
change the shape of the peak at 160 GeV. Similar results are obtained for muons in the 
final state. Detector effects and background processes thus will have little effect on the 
observability of the peak originating from the WW box diagrams. 

A simple method to estimate whether one can hope to observe the peak at 160 GeV in 
Run II is to compare the statistical uncertainty expected for A-fb in a 10 GeV bin centered at 
160 GeV with the variation of the forward-backward asymmetry due to the non-factorizable 
weak corrections in the same region. For 20 fb _1 , the statistical uncertainty in the 160 ± 
5 GeV bin is found to be SApB{stat) ~ 0.016 per lepton channel and experiment, whereas 
the non-factorizable weak corrections change A^b by about 0.003. It will thus be difficult 
to observe the threshold effect in the WW box diagrams in Run II. 

For the definition of cos 8* given in Eq. (0), A^b — for pp collisions. The easiest 
way to obtain a non-zero forward-backward asymmetry at the LHC is to extract the quark 
direction in the initial state from the boost direction of the di-lepton system with respect to 
the beam axis ||47|| . The cosine of the angle between the lepton and the quark in the 
rest frame is then approximated by 



Pz{l + l ) m(l+l-)^m 2 (l+l-) +p 2 T {l+l-) 



p + (r)p-(l + ) -p-(r)p + (l + ) . (40) 



At the LHC, the sea - sea quark flux is much larger than at the Tevatron. As a result, 
the probability, f q , that the quark direction and the boost direction of the di-lepton system 
coincide is significantly smaller than one. The forward-backward asymmetry is therefore 
smaller than at the Tevatron. Events with a large rapidity of the di-lepton system, y(l + l~), 
originate from collisions where at least one of the partons carries a large fraction x of the 
proton momentum. Since valence quarks dominate at high values of x, a cut on the di-lepton 
rapidity increases f q , and thus the asymmetry |47|] and the sensitivity to the effective weak 



mixing angle. In the following we therefore impose a 

|y(/ + r)|>l (41) 

cut in all numerical calculations of the forward-backward asymmetry at the LHC. 

For muons in the final state we impose the pr and pseudo-rapidity cuts listed in Eq. (p0|). 
In the e + e~ case, we allow one of the electrons to be in the range \r)(e)\ < 4.9, whereas the 
other electron is required to be within \i](e)\ < 2.5. This takes into account the possibility of 
using the forward calorimeter in ATLAS for electron identification. The standard rapidity 
coverage of the ATLAS and CMS detectors (see Eq. fl3TD) for leptons, < 2.5, is known to 
significantly reduce A FB @- ^ n addition, it results in a reduction of the total cross section 
in the Z pole region by roughly a factor 5. Combined, these effects greatly reduce the 
chances for a precise measurement of the weak mixing angle at the LHC. As demonstrated 



in Ref. [24|, a large fraction of the sensitivity lost can be recovered if one can make use of 



the forward calorimeter to detect one of the electrons. 
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AA FB at the LHC for 40 GeV < m(7+r) < 200 GeV is shown in Fig. 0. In the 
muon channel, AA FB is smaller than 0.001 in magnitude over the entire mass range consid- 
ered and the peak caused by the threshold effects associated with the WW box diagrams 
is significantly washed out (dashed line). For electrons, on the other hand, the peak at 
m(e + e~) « 160 GeV is quite pronounced and AA FB is roughly a factor 2 larger than in 
the muon case. The dotted line shows the result for electrons, taking into account elec- 
troweak background processes. To reduce the it, W + W~ and ZZ backgrounds we require 
that f T < 20 GeV, and that no jets with px(j) > 50 GeV and \r)(j)\ < 5 are observed. As 
at the Tevatron, the ZZ background is negligible. However, since the it cross section at the 
LHC is more than a factor 100 larger than at the Tevatron, the it background is much more 
important at the LHC and cannot be neglected. For m(e + e~) > 100 GeV, the electroweak 
background processes are seen to significantly modify AAfb, however, without affecting the 
shape of the peak. The expected statistical errors for the forward backward asymmetry in a 
10 GeV bin centered at m{l + l~) = 160 GeV for an integrated luminosity of 100 fb _1 (1 ab _1 ) 
are 5A FB (stat) = 0.0036 (5A FB (stat) = 0.0011) in the electron, and 5A FB (stat) = 0.0062 



(5A FB (stat) = 0.0020) in the muon channel. From Fig. |K] it is then clear that observation 
of the peak caused by threshold effects associated with the WW box diagrams at the LHC 
will also be quite difficult. 

In the Z peak region, the forward-backward asymmetry, A FB , provides a tool to measure 
sin 2 6 l eS . In this region, the forward-backward asymmetry can to a very good approximation 
be parameterized by |4£| 



A FB = b (a - sin 2 6 l cS ) , (42) 

both in the Born approximation and including 0(a) electroweak corrections. The parameter 
b controls the sensitivity of A FB to the effective weak mixing angle. The values of the 
coefficients a and b in the EBA, and the shifts introduced by the QED and the non-universal 
weak corrections are listed in Table |I| for the integrated forward backward asymmetry in 
the region 75 GeV < m(/ + /~) < 105 GeV. The coefficients a and b including the full 0(a) 
electroweak corrections are given by 

a°(° 3 ) = a Q ED + Aa weak , b°^ = 6 QED + A6 wcak . (43) 

Here 



a 



QED = a EBA + Aa QED ; 6 QED = fe EBA + A6 QED ^ 



denote the parameters obtained when only the QED corrections and the factorizable elec- 
troweak corrections in form of the effective Born approximation are taken into account. 



Since previous measurements of A FB at the Tevatron |44| have corrected for detector 
effects, we do not impose any cuts or lepton identification requirements when extracting 
a and b for pp collisions at y/s = 2 TeV. For the LHC we impose the cuts and lepton 
identification requirements described above and in Sec. [Ill A| . In order to obtain the results 
listed in Table [H], we have varied the Higgs boson mass between 75 GeV and 350 GeV, 
corresponding to a variation of sin 2 9 l eS between 0.23149 and 0.23225. 

Table [H] shows that the non-factorizable weak interactions have only a small effect on a 
and b. QED corrections, on the other hand, have a significant impact. In particular they 
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reduce the sensitivity of A F b to the effective weak mixing angle. The rather large differences 
between the coefficients a and b for e + e~(7) and final states at the LHC are due 

to the different rapidity coverage assumed for electrons and muons. The sensitivity of A F & 
to sin 2 9 l cS at the Tevatron is significantly higher than at the LHC [Dj] . 

At the Tevatron in Run II, one expects to measure sin 2 9 eS with a precision of 0.0005 
(0.0006) in the electron (muon) channel, assuming an integrated luminosity of 10 fb _1 ||. 
At the LHC, with 100 fb _1 and the rapidity cuts described above, one hopes to reach an 
accuracy of 0.00014 for the effective weak mixing angle in the e + e~(7) final state p4| . 
Ignoring the 0(a) electroweak radiative corrections would shift sin 2 9 l eS by (2 — 3) 10~ 4 
((3 — 5) 10~ 4 ) towards smaller (larger) values at the Tevatron (LHC). The shift is of the 
same size (Tevatron) or larger (LHC) than the expected experimental uncertainty. It will 
thus be necessary to take 0(a) corrections into account when one extracts the effective weak 
mixing angle from Drell-Yan production in future Tevatron or LHC experiments. If only the 
non-universal weak corrections are neglected, the calculated and the true value of sin 2 6 l cS 
deviate by (2 — 3) 10~ 5 ((1 — 2) 10~ 5 ). Non-universal weak corrections thus have a small 
or negligible effect on the effective weak mixing angle extracted from the forward-backward 
asymmetry. 

Analogously to the invariant mass distribution, one expects that the genuine weak 
corrections to A f -q become large for high di-lepton invariant masses. The forward-backward 
asymmetry at the Tevatron for m(l + l~) > 200 GeV is shown in Fig. |Tl|. The QED correc- 
tions gradually increase in size with increasing invariant masses. For m(l + l~) = 200 GeV 
(m(l+l-) = 600 GeV) they decrease A FB by 0.007 (0.012), i.e. by about 1.2% (2.0%). The 
non-factorizable weak corrections further reduce A F b for the range of masses shown and in- 
crease steadily in size with m(l + l~). For m(l + l~) > 600 GeV they are larger than the QED 
corrections. The difference in shape between the forward-backward asymmetry for electrons 
and muons in the final state at large di-lepton invariant masses is due to the different mo- 
mentum resolution for the two particles. For high energy electrons, the resolution a/E is 
becoming independent of the momentum. For high energy muons, on the other hand, the 
resolution is proportional to the momentum of the particle, and a/p becomes of 0(1) for 
momenta in the several hundred GeV region. The momentum of a high energy muon thus 
is easily mis-measured by a factor 2 or more, thus modifying A F &. 

The forward-backward asymmetry for di-lepton invariant masses between 200 GeV and 
2 TeV at the LHC is shown in Fig. [T2|. For a invariant mass of 2 TeV, the weak cor- 
rections are about a factor two larger than the QED corrections. The full 0(a) electroweak 
corrections reduce A FB by 0.025 (0.04) for electrons (muons), i.e. by about 5% (8%). While 
the forward-backward asymmetry in the high di-lepton invariant mass region at the Tevatron 
is nearly constant, it increases significantly with m(l + l~) at the LHC. For growing di-lepton 
invariant masses, the average fraction x of the proton momentum carried by the quarks 
increases, and, as discussed before, this leads to larger values of A FB . The large difference 
in the magnitude of A FB for electrons and muons is due to the fact that one of the electrons 
may have rapidity up to \t]\ = 4.9, whereas both muons have to be in the range |//(^)| < 2.5. 

While the non-factorizable weak corrections to the forward-backward asymmetry increase 
in magnitude with energy, they are significantly smaller than for the di-lepton invariant mass 
distribution. This is mostly due to an accidental cancellation of the box diagrams and the 
W contributions to the vertex corrections 114 . 
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IV. CONCLUSIONS 



Drell-Yan production in hadronic collisions is an important process. With the anticipated 
large data sets from Run II of the Tevatron and the LHC, it is vital to understand higher 
order QCD and electroweak radiative corrections. In this paper we have presented a calcu- 
lation of the 0(a) corrections to pp — > 7, Z — > l + l~ based on the complete set of one-loop 
Feynman diagrams contributing to di-lepton production. In addition, our calculation takes 
into account the effects of the (9(g 4 m 2 /M^) corrections on sin 2 9 l eS and M w . The calcula- 
tion is based on a combination of analytic and Monte Carlo integration techniques. Lepton 
mass effects are included in the approximation where only mass singular terms originating 
from the collinear singularity associated with final state photon radiation are retained. The 
ultraviolet divergences associated with the virtual corrections are regularized using dimen- 
sional regularization and the on-shell renormalization scheme | 20fl . A previous calculation 
took the 0(a) QED corrections into account JT3[, but ignored the effects of non-universal 



weak corrections. 

The electroweak 0(a) corrections to neutral-current Drell-Yan production naturally de- 
compose into QED and weak contributions which are separately gauge invariant. The QED 
corrections can be further divided into gauge invariant subsets corresponding to initial and 
final-state radiation. The collinear singularities associated with initial state photon radia- 
tion are universal to all orders in perturbation theory and can be absorbed by a redefinition 
of the parton distribution functions. The weak corrections can be written in form of mo- 
mentum dependent effective vector and axial-vector couplings, and contributions from WW 
and ZZ box diagrams. 

Since the phenomenological implications of the QED corrections were discussed in an 
earlier paper ||13|| , we concentrated on the effect of the weak corrections on the di-lepton 
invariant mass distribution and the forward-backward asymmetry. The weak corrections 
were found to enhance the cross section in the Z peak region by about 1%. In contrast, 
QED corrections reduce the cross section in this region by several percent with the exact 
amount depending on the flavor of the charged lepton in the final state, and the lepton 
identification requirements imposed. Since the weak corrections are non-uniform in the 
vicinity of the Z peak, they are expected to shift the Z boson mass extracted from data by 
several MeV. Comparison with the expected statistical precision of about 0.2% (0.05%) for 
the cross section in the Z peak region for 2 fb _1 (10 fb _1 ) at the Tevatron (LHC) shows 
that it will be necessary to take the weak corrections into account when one uses observables 
such as the Z boson cross section or the W/Z cross section ratio to confront data and SM 
predictions, or when calibrating detector components using Z data. 

The non-universal weak corrections were found to have only a small effect on the forward- 
backward asymmetry for m(/ + /~) < 200 GeV. However, threshold effects associated with the 
WW box diagrams lead a characteristic peak in A^b at m(l + l~) « 160 GeV. Unfortunately, 
the size of this footprint of non-factorizable weak corrections is small, and thus will be 
difficult to observe both at the Tevatron and the LHC. 

In the Z peak region, the forward-backward asymmetry provides a tool to measure the 
effective weak mixing angle. Electroweak corrections were found to shift sin 2 9 l cS by an 
amount similar to or larger than the uncertainty expected in future Tevatron and LHC 
experiments, and thus cannot be neglected when extracting the effective weak mixing angle 
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from data. The non-universal weak corrections, however, contribute only 2.5 — 10% to the 
shift. 

The non-factorizable weak corrections to the di-lepton invariant mass distribution and the 
forward-backward asymmetry were found to increase rapidly with m(l + l~). This is due to the 
presence of Sudakov-like electroweak logarithms of the form \n(m(l + l~) / My) (V = W, Z). 
While these corrections are of moderate size (typically a few percent) for di-lepton invariant 
masses accessible at the Tevatron, they become very large for masses in the TeV region 
which play an important role in new physics searches at the LHC For m(l + l~) = 2 TeV, the 
nonfactorizable weak corrections reduce the cross section at the LHC by about 12%. When 
QED corrections are taken into account, the differential cross section may be reduced by 
as much as 40% (see Fig. ||). The strong increase of the non-factorizable weak corrections 
with the di-lepton invariant mass requires that these corrections be resummed. Several 
calculations of the resummed cross section for fermion pair production in e + e~ collisions 
have been carried out recently P5| , however, no such calculation exists yet for di-lepton 
production in hadronic collisions. 

While the electroweak radiative corrections to the di-lepton invariant mass distribution 
at high values of m(l + l~) are very large, they were found to be numerically less important 
in the forward-backward asymmetry. 



ACKNOWLEDGMENTS 

We would like to thank S. Dittmaier, S. Eno, H. Frisch, B. Knuteson, K. Sliwa and 
J. Womersley for stimulating discussions. We also thank T. Hahn for cross checking our 
results for the weak one- loop corrections with FeynArts, FormCalc and LoopTools |50]. One 



of us (U.B.) is grateful to the Fermilab Theory Group, where part of this work was carried 
out, for its generous hospitality. This work has been supported in part by Department of 
Energy contract No. DE-FG02-91ER40685 and NSF grants PHY-9600155 and PHY-9970703. 



20 



APPENDIX A: THE Z DECAY WIDTH 



The total Z decay width Tz is obtained from the sum over the partial decay widths into 
fermion pairs as follows 



v - J2 T ff ■ 



At lowest order in perturbation theory the partial decay widths are given by 



r (o) 



NfT y/l-Anf (l+2/* / )«J + (l-4/i / ) 



with the color factor N9 = 1, 3, / = I, q 



(Al) 



(A2) 



and 



I'f 



m f 
Ml 



(A3) 



The fermionic partial decay widths including electroweak and QCD radiative corrections can 
be expressed in terms of the effective coupling constants gf' V , gf' A and the Z wave function 
renormalization contribution Il z of Eqs. flTSP and flT8|): 



,(0+1) 

// 



NfT 



yi - 4/// 

f^efl z (M; 



1 + 2^) ^• / (M|)- + (l-4 /U/ ) g^ ,T (M 



zj, 



x + 



A, 



c 



'QCD 



(A4) 



The QED corrections, 



°QED 



3aQ 2 f 

47T 



(A5) 



are at most 0.17% of the lowest-order decay width. The QCD corrections for massless 
hadronic final states have been calculated in |5T| , |52| and can be parametrized in the form 
(a s = a s (Mf)) 



'QCD 



+ 1.405 



12.8 



a c 



7T 



Qj aa s 



4 7T 2 



(A6) 



Here, the 0(aa s ) term has also been added although it is not a pure QCD contribution. 

For 6-quarks and r leptons it is important to take into account mass effects in the 
calculation of the electroweak and QCD corrections. From a comparison with a calculation 
for massive external fermions one finds 



r 



(0+1) 
bb 



r^ +1) (/i fe = 0) - 0.0088 GeV , T<? T +1) = T^ T +1 \fi T = 0) - 0.00018 GeV . (A7) 
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APPENDIX B: RENORMALIZED SELF ENERGIES AND FORM FACTORS 



The renormalized self energies fl x (q 2 ) (X = j,Z,jZ) of the neutral vector bosons are 
given by 



^(q 2 )-q 2 W(0), 
S z (g 2 ) - fte£ z (Mf ) + (q 2 - Mf ) 
s w £^(0)\ 



(Bl) 



'<5Mf 



Ml 



'c w Mf 



1 



rr(o) 



S^(g 2 ) = S^(g 2 )-S^(0)-g 2 



5M| 
Mf 



5M 2 



w _ 2 s_ w E^(0) 



c 



,.. M 2 



with n 7 (0) = (9S 7 /9g 2 )|q2 =0 and the mass renormalization constants 



5M Z = Tie £ Z (M|) - 



[E^(M 2 )] 2 



5Mj v = Kei: w (Mj v 



(B2) 
(B3) 

(B4) 



Mf + £T(Af|) / 

where <5Mf is calculated via iteration. E x (g 2 ) (X = 7, Z, jZ, W) denote the unrenormalized 
self energies as the transverse coefficients in the expansion 



9^ X (q 2 ) + ^ ^ X (q 2 )-^(q 



(B5) 



The terms proportional to q^q v yield contributions proportional to m 2 in the ON-SHELL 
amplitudes and hence vanish in the limit m/ — ■> 0. Explicit expressions for the unrenor- 
malized vector boson self energies S x (X =j,Z, jZ, W) and the renormalized form factors 
F^' 7),/ , G% ,T),/ are provided in Appendix B and C.l of Ref. ||. 

Higher-order (irreducible) corrections associated with the p parameter can be taken into 
account by performing the replacement 



in Eqs. (|B2| ) and 



5Mf 
Mj~ 
with 



M 2 W 



5M Z 
Mf" 



6M 2 , 
Ml 



+ Ap 



HO 



Ap HO = 3 G » m2t 



8tt 2 V2 



G ^ Ap^(m 2 /M 2 H ) + Cl + c 2 



8n 2 y/2 



71 



(B6) 



(B7) 



The coefficients C\ and c 2 describe the first and second-order QCD corrections to the leading 
G^m 2 contribution to the p parameter, calculated in |54| and |55]| , respectively Their explicit 
expressions can be found in Ref. |3(J (Eqs. (83,84)). a s (m 2 ) is given by the following relation: 



a s (m 2 ) 



12tt 
~23 



In 



m t 



+ 



12tt 



23a s (M 2 ) 



(B8) 



The function Ap( 2 \m 2 / Mfj) describes the leading two-loop electroweak corrections to the p 



parameter and can be found in Ref. [56 
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APPENDIX C: THE BOX CONTRIBUTION 



Two different topologies of weak box diagrams contribute to the process qq — > 7, Z — > 
which we denote as 'direct' and 'crossed' box diagrams. The corresponding contri- 
butions to the one-loop matrix element are labeled by D and C, accordingly. Since all 
fermion mass effects, except the logarithmically divergent terms associated with the final 
state collinear divergences, are neglected, Higgs boson exchange does not contribute, and 
we only have to consider box diagrams involving Z and W ± exchange. For down (up) type 
quarks, the crossed (direct) box diagram with W ± exchange does not contribute. 

The contribution of the box diagrams to the differential cross section dew of Eq. ( |TTD 
thus can be decomposed in the following way 



da 



, „ 327TCT r . 

dP 2 f —5— He J2 [ B d(s, t, My) + B a (s, u, M v ) 



box 



(CI) 



v=z,w 



with 



B D (s, t, My) = k+ (s + tf [2D° 2 + t{D{ + Dl + D\ + D\ + Df + D 12 )} 

+ k v t 2 [8D° 2 + t (D\ + 2D\ + D\ + 2{Dl + Df + D 12 )) - 2sD 
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2 J J 



Bc(s, U, My) = —B D (§, t, My) With t <-> U, Ky K 



where the four-point functions are denoted by D\ = D{ (£, 0, My, 0, My) and 
1 



'- w 



K 



W 



8si 
0. 



(vi + a t )(v q + a q ) — ^ + — 



(vf + 3 Vl af)(v 3 q + 3v q a 2 g ) + (af + 3a l vf)(a 3 q + 3a q v 2 q ) 

QlQq \/ 2 1 2\t 2 1 2\ 1 a 

—z— + o,)(u 9 + a q ) + Aviv q aiG 
(vf + 3viaf)(vg + 3v q a 2 ) - (af + 3a t vf)(a 3 q + 3a q v 2 ) 

QlQq I"/ 2 1 2\t 2 1 2\ ^ 



|s-M 2 | 2 



|s-M 2 | : 



(C2) 



(C3) 



with the vector and axial vector couplings, Vf,cif(f = l,q), of Eq. (E|) and M c being the 
complex Z boson mass. The explicit decomposition of the vectorial and tensorial four point 
functions 



16vr 2 
d A k 



Df,,f,u(i, mi,M v ,m q , M v 



kfi: k^k v 



(2tt) 4 (k 2 - mj)((k - k + ) 2 - M v )((k + p- k + ) 2 - m 2 )((k + k.) 2 - My) 



(C4) 



with 
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Dp = -k +tl Dl + {p- k+)^D\ + k^D\ , 

= k +tl kj rU D\ + {p — k + )^(p — k + ) u D\ + k_^k_ u D2 + g^D^ 

- [k +fJ ,{p - k+) u + k +u (p - fc + ) /Lt ]£>2 2 - (k +ti k- v + k^k^^D] 3 

- [k+n(p ~ k+) v + k +u (p - k + )^Df + [k^{p - k + ) v + k_ u (p - k+^Df , (C5) 

can be found in Ref. |57]]. The expressions for the crossed box diagram, 
D^ u {u, mi, M v , m g , My) can be obtained by replacing p — > p in Eqs. ( P4| ) and (|C5|). 
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TABLES 



TABLE I. The cross section ratios K EW = /a EBA and K QED = a ^ ED /a EBA for 

pp /+r( 7 ) (J = e, /x) at = 2 TeV with 75 GeV < m(Z+r) < 105 GeV. Shown are the 
predictions with and without taking cuts and lepton identification requirements into account. 





with lepton id. 


without lepton id. 




requirements 


requirements 


channel 


^QED 


J(EW K QED 


pp — ► e + e~(7) 


0.988 0.978 


0.949 0.939 




0.936 0.926 


0.981 0.971 



TABLE II. The coefficients a and b defined in Eq. ( |42| ) in the EBA and the shifts intro- 
duced by QED and weak corrections for the integrated forward-backward asymmetry in the region 
75 GeV < m(l + l~) < 105 GeV at the Tevatron and the LHC. The values listed for the Tevatron 
are obtained without imposing any cuts or lepton identification requirements. At the LHC, the 
cuts discussed in Sec. QIC have been imposed, together with the lepton identification requirements 



listed in Sec. Ill A 



a) Tevatron 



final state 


a EBA 


Aa QED 


Aa weak 


^EBA 


A6 QED 


foweak 


e + e~(7) 


0.24585 


0.00221 


-0.00016 


3.408 


-0.408 


0.026 


M + M _ (7) 


0.24585 


0.00094 


-0.00001 


3.408 


-0.171 











b) LHC 








final state 


a EBA 


Aa QED 


Aa weak 


^EBA 


A6 QED 


A ^weak 


e + e~(7) 


0.24797 


0.00284 


-0.00020 


1.618 


-0.276 


0.013 




0.25072 


-0.00012 


0.00037 


0.724 


-0.050 


0.013 
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FIGURES 





FIG. 1. Born and higher-order weak contributions to qq — ► 7, Z — > Z + Z in symbolic notation. 
The dark blob indicates the inclusion of all 1PI contributions to the photon and Z propagators. 




FIG. 2. Box diagrams contributing to qq — > Z + Z 
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FIG. 3. The relative corrections to the total cross sections for uu — > e + e and dd — > e + e a) 
in the vicinity of the Z resonance, and b) at high parton center of mass energies. 



30 



0.04 



0.02 



m 

O El, 

<! 



m 
< 



0.00 - 



-0.02 — 



-0.04 



-0.06 



n — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — r 



a) solid: uu^e + e , d = d 



(o) - 



dots: dd^e + e , a — a^ 
dash: uu->e + e~, c? = cr EBA 



dashdot: dd->e + e , a = d 



EBA 




50 60 70 80 90 
VI (GeV) 



100 



0.06 



0.04 



0.02 



0.00 



-0.02 



-0.04 



" b) solid: uu -> e + e 


1 1 1 


= <7<°> ^ 




u 


= a (0) '. 


dash: uu^e + e", 


^0 = 


CT EBA 


dashdot: dd^e + e~ 




= ° EBA " 




















- 1.' \ ~ ' 














, , 1 





200 400 600 800 1000 
VI (GeV) 



FIG. 4. The relative corrections to the forward-backward asymmetry at parton level for 
uu — > e + e _ and dc? — > e + e~ a) in the vicinity of the Z resonance and b) at high parton cen- 
ter of mass energies. 
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FIG. 5. The relative corrections to the forward-backward asymmetry at the parton level in the 
region around the VF-pair production threshold, y/§ = 2Mw, for a) uu — > e + e~ and b) dd — > e + e~~ . 
The solid line shows ^4fb — ^fb"^ wrien the full set of Feynman diagrams contributing to the 
non-photonic weak corrections is taken into account. The dashed (dotted) lines show ^4fb — ^4fb A 
when the W (Z) box diagrams are disregarded. 



32 



1.01 



1.00 



0.99 



0.98 



i i r 



j i i_ 



50 



i i i r 

pp ^i + r(7) 

Vs = 2 TeV 

solid: l = e 
dash: l = /x 




J I I L 



100 150 
m(l + r) (GeV) 



200 



FIG. 6. The ratio [da°^ / dm{l + l~ )]/[da^ EB /dm(l+l- )] function of the di-lepton invari- 
ant mass at the Tevatron in the Z peak region. a°^ a > denotes the full NLO cross section, and 
^.QED represents the cross section which includes the factorizable electroweak corrections in form 
of the effective Born approximation together with the 0(a) QED corrections. The cuts and lepton 



identification requirements imposed are described in Sec. Ill A 
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FIG. 7. The ratio [da/dm(l + l~)]/[da EBA /dm(l + l~)] as a function of the di-lepton invariant 
mass for a) pp — * e + e~( , y) and b) pp — > /i + \i~ (7) at = 2 TeV. The dashed lines show the ratio 
of the complete C(a 3 ) electroweak and the EBA differential cross section. The solid lines display 
the corresponding ratio for the case where only the 0(a) QED corrections and the factorizable 
electroweak corrections in form of the effective Born approximation are taken into account. The 
cuts and lepton identification requirements imposed are described in Sec. [II A. 
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FIG. 8. The ratio [da/dm(l + l~)]/[da EBA /dm{l + l~)] as a function of the di-lepton invariant 
mass for a) pp — * e + e~(j) and b) pp — * /i + ^~(7) at = 14 TeV. The dashed lines show the ratio 
of the complete C(a 3 ) electroweak and the EBA differential cross section. The solid lines display 
the corresponding ratio for the case where only the 0(a) QED corrections and the factorizable 
electroweak corrections in form of the effective Born approximation are taken into account. The 
cuts and lepton identification requirements imposed are described in Sec. [II A. 
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FIG. 9. The difference A FB (full EWK) - A FB (QED) for pp -> at Vs = 2 TeV. 

Ape (full EWK) denotes the forward-backward asymmetry calculated taking the full 0(a) elec- 
troweak corrections and the 0(g 4 m^/My V ) corrections into account. Afb(QED) only includes 
the 0(a) QED corrections, in addition to the factorizable corrections absorbed in the EBA. The 
solid and dashed lines show the results for electron and muon final states, respectively. The dot- 
ted line shows the difference in the asymmetry taking the pp — > W + W~ —* e + e~p' T background 
with f T < 20 GeV in ^FB(full EWK) into account. Additional cuts and the lepton identification 



requirements imposed are described in Sec. Ill A 
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FIG. 10. The difference A FB {Ml EWK) - A FB (QED) for pp -» Z+r( 7 ) at y/s = 14 TeV. 
4.pB(full EWK) denotes the forward-backward asymmetry calculated taking the full 0(a) elec- 
troweak corrections and the 0(g 4 m^/My V ) corrections into account. Afb(QED) only includes 
the 0(a) QED corrections, in addition to the factorizable corrections absorbed in the EBA. The 
solid and dashed lines show the results for electron and muon final states, respectively. The 
dotted line shows the difference in the asymmetry taking the pp — > W + W~ — > e + e~^> T and 
pp — > it — > e + e~f T bb background in .ApB(full EWK) into account. The cuts and lepton identifica- 
tion requirements imposed are described in Sees. [Ill A| and IIIG 
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FIG. 11. The forward-backward asymmetry as a function of m{l + l~) for a) pp — > e + e~{^) and 
b) pp — * fi + /j,~(^f) at y/s = 2 TeV. Shown are the asymmetry in the EBA (solid lines), including 
pure QED corrections in addition to those corrections which are part of the EBA (dashed lines), 
and the asymmetry taking the complete set of 0(a) electroweak corrections and the 0(g 4 m^ /Afjjp) 
corrections into account (dotted lines). The cuts and lepton identification requirements imposed 
are described in Sec. Ill A| . 
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FIG. 12. The forward-backward asymmetry as a function of m{l + l~) for a) pp — > e + e~{^) and 
b) pp — > /i + //~(7) at y/s = 14 TeV. Shown are the asymmetry in the EBA (solid lines), including 
pure QED corrections in addition to those corrections which are part of the EBA (dashed lines), 
and the asymmetry taking the complete set of 0(a) electroweak corrections and the 0(g 4 mf /Afjjp) 
corrections into account (dotted lines). The cuts and lepton identification requirements imposed 
are described in Sec. Ill A| and [III C] , 



39 



